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Identification of the ionizing source of NGC 2024* 
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Abstract. We propose the late-O, early-B star IRS2b as the ionizing source of the Flame Nebula (NGC 2024). It 
has been clear that such a hot, massive star must be present in this heavily obscured region, and now it has been 
identified. New near-infrared photometry shows that IRS 2b is the most luminous and hottest star in the young 
star cluster embedded in the center of NGC 2024. The near-infrared observations (5' x 5') cover ~90 % of the 
H n region detected in radio continuum radiation, making the probability very low that the ionizing star is not 
present in the field. A K-band spectrum of IRS2b obtained with ISAAC on the Very Large Telescope indicates 
that the spectral type of IRS2b is in the range 08 V - B2 V. Additional arguments indicate that its spectral type 
is likely closer to 08 than to B2. The corresponding amount of ionizing radiation is consistent with published 
radio continuum and recombination line observations. 
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O , 1. Introduction 



Just east of ( Ori and north-east of the Horsehead Nebula, 
NGC 2024 (Flame Nebula) appears in visible light as a 
bright nebula of which the central part is obscured by a 
thick dust lane. The central region of NGC 2024 is also 
a bright sourc e of radio continuum emission and rec om- 
bination lines l|Krugel et al.lll98"3lBarnes et al.l fl989). in- 
dicating the presence of an ionizing star of spectral type 
09-09.5. The heavy obscuration b y dust is the reason 
why, contrary to M42 l)O'dellll200lj) and many other vis- 
ible H II regions, the exciting star(s) of NGC 2024 have 
not been identified. Near- and mid-infrared observations 
have revealed the p resence of a young star cluster in the 
core of NGC 2024 fearnes et al.lll98a lLada et alJll99ll 
IComeron et al. I ll996t lHaisch et alJ \200% ~ Several candi- 
dates f or the ionizing source of NGC 2024 h ave been pro- 
posed l|Grasdalenlll974l iBarnes et alJll989|h but none of 
these candidates are able to produce the observed radio 
continuum emission. 

We are carrying out a systematic survey of the stellar 
content of compact and ultra-compact H n regions, with 
the aim to detect and study deeply embedded, newly-born 
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* Based on observations collected at the European Southern 
Observatory at La Silla and Paranal, Chile (ESO programmes 
62.H-0443 and 64.H-0425) 



massive stars. Our ultimate goal is to better understand 
the earliest phases in the life of the most massive stars. In 
this context we have obtained deep near-infrared images 
of a sample of com pact H n regions, including NGC 2024 
l|Kaper et al.lbOO.'j) . Subsequent K-band spectroscopy of 
the candidate ionizing stars h as resulted in t he determi- 
nation of their spectral types l)Bik et alJl2003lL 

Here we describe how we have identified the long- 
sought-for ionizing star of NGC 2024. In Sect.|2]we present 
near-infrared, narrow-band images of this region, as well 
as a K-band spectrum of the candidate ionizing star 
IRS2b. In Sect.|31we derive the physical properties of this 
star. In the last section we present our conclusions and 
compare the derived properties of IRS2b with alternative, 
previously proposed candidates. 



2. Observations 

Near-infrared images of the central region of NGC 2024 
were obtained with SOFI mounted on ESO's New 
Technology Telescope (NTT) at La Silla on February 6, 
1999; the seeing was 0.7". We used narrow-band filters 
centered on strong nebular emission lines (P/3 1.28 /an, 
H2 2.12 /im, Br7 2.16 /an), and two narrow-band contin- 
uum filters in the J and K-band (1.21 and 2.09 /im). The 
latter were chosen to measure the J and K magnitude of 
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the embedded stars, avoiding the contamination by nebu- 
lar emission lines. 

Nine frames of 20 seconds each were taken on source 
and 9 frames with a 2 second exposure per frame were 
obtained to measure the sky background. The observa- 
tions were reduced using standard infrared imaging data 
reduction procedures with IRAF. 

Pointsources were detected by adding the frames in 
the J and K c ontinuum filters, and running DAOPHOT 
l|Stetsonlll987h on the resulting frame. The photometry 
was obtained as follows: aperture photometry with a large 
aperture was performed on the image of the standard star 
and of bright, isolated stars in the image field. This allowed 
us to set up a network of bright secondary standards in the 
image field. Then, aperture photometry with a small aper- 
ture (3 pixel in radius), adequate for our rather crowded 
field, was performed on all the stars in the field. The mag- 
nitudes of these stars were determined taking those of the 
secondary standards as a reference. 

We note that for red sources ((J-K) > 2 mag) the 
narrow-band J and K magnitudes differ significantly from 
the broad-band J and Ks magnitudes. This difference has 
been estimated by multiplying various energy distribu- 
tions with the response curves of the filters, representing 
a range in (J-K) from to 7. The typical errors for (J- 
K) « 5 on this correction are 0.2 and 0.07 mag in J and 
K, respectively. The J and K-band magnitudes used in 
this paper are based on narrow-band observations. Only 
for the construction of the Hertzsprung-Russell diagram 
(HRD, Fig. 2J) the narrow-band magnitudes are converted 
into broad-band magnitudes. 

Medium-resolution (R ~ 8, 000) K-band spectra of 
some stars in this region were taken with ISAAC and 
ESO's Very Large Telescope (VLT) at Paranal, Chile on 
March 20, 2000. The spectra were dark- and flat-field cor- 
rected, and wavelength calibrated using standard reduc- 
tion procedures. In order to correct for the sky background 
the object was "nodded" between two positions on the slit 
(A and B) such that the background emission registered at 
position B (when the source is at position A) is subtracted 
from the source plus sky background observations at posi- 
tion B in the next frame, and vice versa. Telluric absorp- 
tion lines were removed using the telluric standard star 
HD 39908 with spectral type A2 V observed under identi- 
cal sky conditions. The only photosphcric line (Br7) in the 
spectrum of the telluric star needs to be devidcd out first. 
It turns out that the best result is achieved when first the 
telluric features are removed from the K-band spectrum 
of the telluric standard using a high resolution telluric 
spectrum (obtained by NSO/Kitt Peak). This spectrum 
is taken under completely different sky conditions, so a 
lot of remnants are still visible in the corrected standard 
star spectrum; without this "first-order" telluric correc- 
tion, a proper fit of Br7 cannot be obtained. The Br7 line 
is fitted by a combination of two exponential functions. 
The resulting error on the Br7 equivalent width (EW) of 
our target star is about 5%. 



An L-band spectrum of IRS2b, with a spectral res- 
olution R ~ 1200 was obtained with VLT/ISAAC on 
February 22, 2002. The data reduction was performed 
using an A6 II star (HD 73634) as telluric standard. The 
quality of the telluric standard was not sufficient to obtain 
a quantitative measurement of the hydrogen absorption 
lines in the L band spectrum. 

3. The nature of the ionizing source in NGC 2024 

Fig. ^ shows a composite near-infrared image of the ob- 
scured central region of NGC 2024 (5' x 5'), centered 
on the bright (saturated, K ~ 5) infrared source IRS2. 
Exposures in three narrow-band filters are combined: Bi'7, 
H2 and P/3. The edge of the molecular cloud to the south 
of NGC 2024 is marked by emission produced by excited 
molecular hydrogen: note the clumpy filamentary struc- 
ture. The recombination lines of hydrogen (e.g. Br7) are 
strongest towards the core of the embedded stellar cluster. 
The optical extent of the Flame Nebula is roughly 15' x 15', 
but at radio wavelengths NGC 2024 is only slightly larger 
(6' x 6') than the field covered by our near-infrared obser- 
vations. The radio continuum radiation is strongly con- 
centrated towards the central part of NGC 2024, and is 
expected to outline the extent of the region where hydro- 
gen is fully ionized. 

In the following we demonstrate that IRS2b is the best 
candidate ionizing star present in the central star cluster 
of NGC 2024. A K-band spectrum of IRS2b is used to 
confirm its early spectral type. The resulting effective tem- 
perature and luminosity of IRS2b are consistent with the 
amount of observed radio continuum and line emission. 

3.1. The young stellar cluster embedded in NGC 2024 

Due to the severe obscuration by the dust lane in front 
of NGC 2024, its stellar population becomes apparent 
at near-infrared wavelengths. The strong infrared sources 
IRS1-5 are the brightest objects in Fig.f^ Another bright 
source, IRS2b is located just 5" north-west of IRS2 (see 
inset Fig. HI ) The existence of IRS2b was first reported b y 



I.Tiang et all 1 1984h and confirmed by iNisini et" 

ID dial- 

We detect 106 and 210 point sources down to the de- 
tection limit of 17.7 mag in J and of 16.1 mag in K, re- 
spe ctively. Because of the proximity of NGC 2024 (360 
nc. lBrown et alJ|l994h . its small spatial extent, its angu- 
lar separation from the galactic plane, and the obscuration 
of background sources by the molecular cloud, we expect 
that nearly all of the infrared sources present in our near- 
infrared image are p hysically related to the cluster (cf. 
IComeron et "afl ll 9961) . 

The brightest infrared source in the center of the 
NGC 2024 cluster is IRS2, a source that has been held 
responsibl e for the ionizat ion of the H II region since its 
discovery ljGrasdalenlll974|) . However, the physical nature 
of IRS2 is unclear: its K-band magnitude is too bright 
for a main sequence star, indicating a strong infrared ex- 
cess. This excess is likely due to the presence of a dense 
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Fig. 1. Composite of three narrow-band images of the central region of NGC 2024 obtained with NTT/SOFI (blue: 
P/3; red: Br7; green: H2). North is up and East to the left; the image size is 5' x 5' (at the distance of NGC2024 this 
corresponds to 0.5 x 0.5 pc). The upper left panel zooms in on the region around IRS2: about 5" to the north-west 
we find IRS2b, the ionizing star of NGC 2024. Also some previously proposed candidate ionizing stars are indicated. 
The vertical stripes are instrumental artifacts caused by the brightest stars. 



circu mstellar disk; the cen tral star is probably a B-type 
star ijLenorzer et al.ll2003l and below). Can we identify 
another, potentially hotter star that could be responsible 
for the ionization of the Flame Nebula? 

The K vs. (J-K) color-magnitude diagram (CMD) of 
the stars detected in our images is presented in Fig. [5J 
only the brightest objects are shown. T he dotted line ind i- 
cates the position of the main sequence ijAller et al.ll 9 821 , 
for different amounts of visual extinction Ay. A distance 
modulus of 7.8 ±0.45 mag is adop ted, which corresponds 
to 363 ± 75 pc ijBrown et alJll994(i . An anoma lous extinc - 
tion law with Ry = 5.5 is used, as measured bvlLed 1U9 68). 
We used the parameterization of ICardelli et alJ i|l989h to 
describe the shape of the extinction law. This extinction 
law is derived for stars with E(B-V) < 1 and could be 
different for higher amounts of reddening. However, the 
dereddening of OB stars in the near infrared is not very 
sensitive to the precise value of Ry, as Ak = 0.108 x Ay 
for Ry = 3.1 and A K = 0.125 x Ay for Ry = 5.5. 
The slope of the near-infrared extinction law changes very 
slowly with Ry and the intrinsic (J-K) color of OB stars 
is almost constan t with spectral type and luminosity class 
(|Koornneedll983|) . 

Fig.Hshows that the star IRS2b (K=7.57±0.07), when 
dereddened to the main sequence, is the best candidate 
ionizing star in the field (neglecting IRS2) at a position 
consistent with a late O main sequence star. The position 
of IRS2b in the CMD implies Ay = 24.0 ±0.5 mag (Ay = 
28.5 ± 0.5 for Ry = 3.1) if it is a main sequence star. 

3.2. K-band spectral type of IRS2b 

The spectral type of IRS 2b can be determined from its K- 
band spectrum (Fig. Hanson et al.l l|l996l) define five 
K-band spectral classes for O and early-B stars based 
on low-resolution K-band spectra. Two K-band spectral 
classes correspond to the hottest O stars (k03-04 and 
k05-06, with k denoting that the classification is based on 
the K-band spectrum) which show lines of N ill (2.115 /mi) 
and Civ (2.079 /im) in emission and the Br7 (2.166 /im) 
line in absorption. In the third class (k07-08) the ClV 
line is absent and the Hei (2.1128, 21137 /im) lines ap- 
pear in absorption. The fourth spectral class k09-Bl is 
defined as having He 1 and Br7 both in absorption, with 
the equivalent width of Br7 less than about 4 A. This class 
is equivalent to a Morgan-Keenan (MK) spectral type be- 
tween 08 V and Bl V. The fifth class (kB2-B3) shows 
strong Br7 and Hei absorption lines and is equivalent to 
MK type Bl V to B2 V. The line equivalent widths (EW) 
can in principle be used to determine the K-band spec- 
tral type. There is, however, a substantial scatter in the 
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Fig. 2. Color-magnitude diagram of the brightest stars in 
the central region of NGC 2024. The vertical (dotted) lines 
indicate the position of the main sequence for different 
amounts of visual extinction (Ay). The magnitudes are 
obtained through narrow-band observations. The diago- 
nal (dashed) line is the dereddening line of IRS2b, demon- 
strating that IRS2b is the intrinsically brightest star of the 
embedded population. IRS2 (saturated in K) probably is 
a B star with a strong infrared excess due to a circumstel- 
lar disk l|Lenorzer et alJl2003|) . The dashed-dotted line is 
the detection limit in (J-K) following from the detection 
limits in J and K. 



observed EW of Bi'7 as a function of spectral type in the 
calibration stars used bv lHanson et al 

IRS2b has both Hei 2.113 /im and Br7 in absorption, 
with an equivalent width of 0.4 ± 0.1 A and 4.9 ± 0.6 A, 
respectively. The Br7 EW indicates that the K-band spec- 
tral type of IRS2b is in the range k09-Bl to kB2-B3 (MK 
types 08-B2). The Hei line appears in absorption around 
spectral type 07.5, increasing its EW towards later spec- 
tral types. The Hei line decreases in strength again to- 
wards B2/B3. The relatively weak Hei line in the spec- 
trum of IRS2b suggests that its spectral type is cither 
close to 08, or more towards B2/B3. 

An L-band spectrum provides a better diagnostic 
of its spectral type, since it includes hydrogen Pfund 
lines, which are more temperature sensitive than Br7 
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Fig. 3. K-band spectrum of IRS2b, obtained with 
VLT/ISAAC, showing Hei 2.1127, 2.1138 /^m and Br 7 
2.166 /jm in absorption. 

Lenorzer et al.ll2"002j) . The L-band spectrum we obtained 
from IRS2b is of insufficient quality to measure the 
strength of the absorption lines, though it is clear that 
they are present. Although we are not able to provide an 
independent estimate of its spectral type, the L-band spec- 
trum rules out the possibility that IRS2b has an infrared 
excess. Such an excess would be even more dominant in 
the L-band, hindering the detection of photospheric lines. 

3.3. The position of IRS2b in the HRD 

Taking our measurement errors into account, we arrive 
at a (conservative) estimate of the MK spectral type of 
IRS2b in the range between 08 V and B2 V. This corre- 
sponds to a range in T P ff between 34 00 and 22,000 K 
l|Martins et alJ 120021: I Alter et all Il982f) . Note, however, 
that the effective temperature calibration of OB stars is 
still a matter of debate. The derived limits on T e s are not 
very hard. 

For every T e g, the dereddened K-band magnitude can 
be converted into a luminosity (adopting a distance of 
363 pc). For this calculation we use t he relation between 
the bolometric correction and T e g from lVacca et al 1 lll996l) 
and the (V-K) vs. T c g relation given bv lK^ornneeTl|l983|) . 
The uncertainty in the luminosity of IRS2b is due to the 
uncertainty in the distance towards NGC 2024 and to the 
conversion from narrow-band to broad-band magnitudes. 
The observational constraints on the location of IRS2b in 
the HRD are visualized in Fig. 21 (left panel, shaded area). 

Are these constraints on L and T c g of IRS2b consis- 
tent with the observed nebular emission? The radio con- 
tinuum flux of the H II region can be used to estimate 
the ultraviolet flux of the ionizing source of NGC 2024. 
The total radio continuum flux at 1667 MHz derived f rom 
high-resolution radio observations ijBarnes et alJfl989T) is 
63 ± 4 Jy. The number of Lyman continuum photons is 
derived from the observed radio continuum flux under the 
assumption that the H n region is ionization bounded and 
neglecting absorption of UV radiation by dust, implying 
that the number of recombinations is equal to the number 



of ionizations. We find that the number of Lyman contin- 
uum photons is 7. 3 ±1.2 x 10 47 The given accuracies refle ct 
th e error in the distanc e quoted bv lBrown et al 1 ill 9941) . 

iKriigel et alJ l)l982l) performed observations of the ra- 
dio recombination lines H76a and Hc76a. The corre- 
sponding ionized helium fraction is 3 %, which leads to 
2.2 ± 0.1 x 10 46 helium continuum photons. 

The observed number of photons capable of ionizing 
hydrogen (Qo) and helium (Qi) can be compared to pre- 
dictions based on s tellar atmosphere models of OB stars 
l|Smith et alJl2002fl . If we assume that a single hot star 
ionizes NGC 2024, Qo and Qi constrain the location of 
the ionizing star in the HRD. In Fig. 0] (right panel), the 
regions which are able to produce the required amount of 
Qo and Qi arc indicated by vertical and horizontal stripes, 
respectively. These regions overlap with the shaded region 
in the HRD defined by the spectral classification of IRS2b. 

4. Discussion 

As shown in the previous section, an 08 star with charac- 
teristics compatible with those observed in IRS2b should 
be capable of producing the required amount of ionizing 
radiation to explain the degree of ionization of NGC 2024. 
To meet the constraint set by Qi, its spectral type is more 
likely to be late O than early B, which is supported by its 
position in the CMD. Note, however, that the predicted 
value for Qi is not well known and depends on the selected 
stellar model. Although other candidate ionizing stars ex- 
ist, they are not needed to provide a self-consistent solu- 
tion for the ionization of NGC 2024. In the following we 
will discuss some of these alternatives. 

If we correct for foreground extinction (cf. Fig. |2J), 
IRS2b and IRS2 are the two brightest sources. The posi- 
tion of IRS2 in the CMD indicates a st rong infrared excess 
((Lenorzer et al.ll2003llGrasdalerll9'7l . IRS2 is associated 
with the ultra-compact radio source G206. 543-16. 347, but 
its physical nature is unclear. Current models suggest that 
ISR2 is a B star surrounded by a dense circumstellar disk 
responsible for the production of the infr ared exces s. 

IRS1 (estimated spectral type B0.5V. lGarrisonlll968() 
and IR S2 have been cand idate ionizing sources for a long 
time. iBarnes et alJ dl989|) discovered 29 additional near- 
infrared sources in the Flame Nebula. In Figs. and 
HI t he brightest of these infrared sources discovered by 
IBarnes et alJ l|1989l) are plotted, as well as IRS1 and IRS2. 
They propose that IRS1, IRS4 and IRS5, together with 
IRS3, provide a significant contribution to the ionizing ra- 
diation. IRS3 is not shown in the CMD, because it is not 
a single source, but consists of mult iple stars which were 
not resolved bv IBarnes et al.1 l(l989h . As Qi is decreasing 
very rapidly with later spectral type the contribution to 
Qi by the other bright stars (all later spectraltype than 
B0.5V) is negligible. 

The radio continuum emission, however, can only be 
transformed into the number of ionizing photons assuming 
that there is no dust included in the Hn region, and that 
the H II region is ionization bounded. If these assumptions 
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Fig. 4. Le/t;The upper part of the theoretical HRD. The shaded region indicates the position of IRS2b for which 
all observational constraints on L and T c ff are fulfilled. The diagonal area represents the constraint provided by the 
observed K-band magnitude. The K-b and spectrum implie s that T e g < 34000 K. Also shown are the theoretical 
evolutionary tracks and isochrones from lLeieune k, Schaererl l|200ll) . Right: The same region of the HRD as shown in 
the left panel, but now including the regions in the HRD capable to produce Qo and Qi for comparison. 



are not valid, the derived values of Qo and Qi are lower 
limits, so that a star hotter than IRS2b might be required. 
In principle, such a star could be significantly more red- 
dened (Ay > 35 mag, based on our detection limit in 
J) and thus have remained undetected, especially when 
the amount of extinction strongly varies with position. 
However, the visual extinction we measure for the stars 
in the field is in the range between 5.5 < Ay < 25 mag. 
We judge that the probability of selectively obscuring this 
potentially hotter star is low. 

In F ig. |H e volu t ionary tracks and isochrones from 
iLeieune k. Schaererl l|200ll) are plotted and the zero-age 
main sequence (ZAMS) is indicated. The location of IRS2b 
suggests that its mass is between 15 and 25 Mq. Although 
it is not possible to determine the age of a single star, 
to position of IRS2b in 0] is, given the uncertainty in 
the effective te mperature calibration consistent with the 
age proposed b y Blaauw (1991); iBrown et al. I Jl994h and 
IComeron et alJ l|l99fij) . IComeron et al.l l|l99fil ) derive an 



age of 2 x 10 6 years based on the frequency of IR excess in 
very low mass objects in NGC 2024 as compared to that 
in Rho Ophiuchi. 
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